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Recent results from Au+Au collisions at BNL-RHIC en- 
ergy hint at explosive hadron production at the QCD transi- 
tion rather than soft hydrodynamic evolution. We speculate 
that this is due to a rapid variation of the effective potential 
for QCD close to Tc- Performing real-time lattice simulations 
of an effective theory we show that the fast evolution of the 
potential leads to "explosive" spinodal decomposition rather 
than bubble nucleation or critical slowing down. 



It is hoped that heavy ion coUision experiments at 
BNL-RHIC and later at CERN-LHC will provide ev- 
idence for the deconfined phase of QCD Lattice 
gauge simulations predict a weakly first-order deconfine- 
ment transition at a critical temperatm'e of Tc ^ O.GS^/ct, 
where a ^ 1 GeV/fm is the string tension at T = 0. For 
QCD with quarks, the transition could remain weakly 
first-order or become a smooth crossover. Even so, lattice 
data for two or three light flavors |^ indicate a steep rise 
in energy density and pressure with temperature and thus 
a rapid shift of the global minimum of the effective ther- 
modynamic potential for QCD within a narrow temper- 
ature interval centered at Tc. To date, it has been diffi- 
cult to identify unambiguous observables providing com- 
pelling evidence for a phase transition in ultrarelativistic 
heavy ion collisions (URHIC). However, rapid changes in 
the QCD effective potential near Tc and/or supercooling 
followed by spinodal decomposition can have distinctive 
signatures. Their experimental observation would proba- 
bly represent the most direct information which URHIC 
events can provide regarding the QCD phase transition. 
The purpose of this Letter is to consider the likelihood 
of explosive "spinodal" decomposition in RHIC experi- 
ments. 

In a first-order phase transition, the thermodynamic 
potential for T slightly less than Tc exhibits a metastable 
minimum corresponding to the deconfined phase. At 
lower temperatures, this minimum becomes a point of 
inflection, i.e. a "spinodal instability" . In a slowly ex- 
panding system, such as the early universe, the phase 
transition proceeds through the nucleation of bubbles of 
the "true vacuum" state via thermal activation for all 
T < Tc In ultrarelativistic collisions of heavy ions 
or hadrons, expansion is very fast and the metastable 
deconfined state may reach the inflection point and un- 
dergo spinodal decomposition before significant nucle- 



ation has taken place. Calculations using an effective 
chiral model jij suggest that spinodal decomposition is 
the favored mechanism in URHIC. The condensate field 
in each causally connected region "rolls down" to its min- 
imum with large attendant fluctuations in the chiral or- 
der parameter and related hadronic observables. 

Effects similar to those of spinodal decomposition can 
also be expected, independent of the order of the phase 
transition, for any rapidly varying thermodynamic poten- 
tial. In this case, flelds are unable to follow the rapid evo- 
lution of the equilibrium potential, and the system decays 
into any of a large number of unstable pion and sigma 
modes. (By contrast, the true spinodal decomposition 
of a slowly varying potential proceeds via the "critical" 
growth of long-wavelength modes of the order parame- 
ter ||. Here, we shall use the term "spinodal" somewhat 
loosely to describe any rapid decomposition.) This was 
first suggested in |6| in the context of chiral symmetry 
breaking. It was argued that, if expansion is sufficiently 
fast, the system can cool significantly while the chiral or- 
der parameter remains near its restored symmetry value. 
Violent decomposition would follow as a natural conse- 
quence of such a "quench" . Subsequent numerical sim- 
ulations showed that, if chiral symmetry breaking is 
considered independent of confinement, the chiral effec- 
tive potential changes too slowly with temperature for 
such effects to be seen in URHIC. If, however, the effec- 
tive potential changes rapidly with temperature the 
quenched scenario and effective spinodal decomposition 
could be an interesting possibility even for moderate ex- 
pansion rates. 

RHIC {y/s — 130A GeV) has already produced a 
wealth of new results "Dynamical" fluctuations in 
the average transverse momentum, pt, are larger than at 
the SPS {y/s = 17 A GeV). Observed Hanbury-Brown- 
Twiss (HBT) radii from pion interferometry are small, 
~ 5 — 7 fm, and average transverse boost velocities are 
large, ^ 0.6c. In particular, the fact that the HBT 
outward and sideward radii, Ro and Rs, are roughly 
equal for pion pair transverse momenta in the range 
100 < Pt < 500 MeV may indicate explosive decom- 
position. These phenomena all hint at a rapid out-of- 
equilibrium transition. By contrast, if a flrst-order tran- 
sition proceeds at equilibrium, a mixed phase is created, 
and the system should spend a long time at and 
slightly below ijlll Tc- As a result, the outward HBT ra- 
dius Ro should be larger than the sideward radius Rg by 
as much as 50% - 100% 
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We wish to consider the possibihty that the evolution 
of central RHIC events is determined by rapid changes 
in the thermodynamic potential. This would be the case, 
for example, if the chiral phase transition were driven by 
the confinement transition. The familiar rapid change 
of the energy density with temperature would then be 
reflected in equally rapid changes in the thermodynamic 
potential for the chiral fields. To illustrate this point, 
we will perform real time simulations of a weakly first- 
order deconfinement transition using a recently suggested 
effective theory, valid near Tc, in which the free energy 
(or pressure) is dominated by a condensate of Polyakov 
loops i ijl^. The potential for £ is taken as ||l^ 

vw = - + in') + b,T\ (1) 

which is invariant under global Z(S) transformations. 
The form of V(^) is dictated by usual symmetry princi- 
ples The center symmetry of SU{3), Z{3), is broken 
at high T where the global minimum of V{£), £o{T) 1, 
and restored at low T where io{T) 0. 

The assumption that quarks are not important for un- 
derstanding the form of the effective potential for QCD 
is motivated by the results of lattice simulations. Specifi- 
cally, plots of P/Pidoai versus T/Tc are remarkably insen- 
sitive to the number of flavors, TV/ = 0, 2, 3 0. Thus, 
we neglect terms linear in i, which must exist in QCD 
with quarks. Such terms would change the first-order 
transition (for Nf = 0) into a crossover. Nevertheless, 
the fact that the free energy is small below means 
that their numerical effects must be small. We empha- 
size that the results (for rapid expansion) to be described 
below are not driven by the order of the phase transition 
and would not be changed in any significant way by the 
inclusion of such linear terms. 

The coefficients 62, 63, and 64 can be chosen to re- 
produce lattice data for the pressure and energy den- 
sity of pure glue theory for T > Tc In practice, 
to reproduce the pressure and energy density we may 
chose constant 63 « 0.9 and 64 « 15.0. For Nf ~ 3, 
we rescale 64 by a factor 47.5/16, corresponding to the 
increase in the number of degrees of freedom. In the 
spirit of mean field theory, the same constant values of 
63 and 64 are to be employed for T < Tc- In this domain, 
&2(r) w -OMa^{T)/T^ can be set using the finite-T 
string tension cr(T) from the lattice 

The potential in eq. (Q) changes extremely rapidly near 
Tc, see Fig. 1 of ||^. For temperatures even ±2% away 
from Tc, the potential has only a single minimum, which 
corresponds to the deconfined (confined) state. This 
rapid change in the effective potential with temperature, 
which is directly related to the rapid increase of o'(T) 
below Tc seen in lattice data makes "explosive" de- 
composition inevitable. (The rapid variation of b2(T) 
about Tc can also be understood qualitatively by consid- 
ering the above mean field theory for the second-order 



confinement transition, 63 = 0, for the quenched theory 
with two colors. Then, 62(T) = c{T/Tc - 1). The rapid 
rise of the pressure relative to that for an ideal gas about 
Tc P(T)/Pid(T) = 5^(T), implies that c is large.) 

To complete the effective theory, we add a kinetic term 
for £ and allow coupling to a chiral field, (f>, describing the 
usual pions and the sigma meson ||^. Thus, 

N 

c ^ + -^\d,efT' - v{e) - —<t>^ \efT^ . (2) 
9 ^ 

Crj, is the standard Lagrangian for the chiral fields, see 
e.g. The coefficient of the kinetic term for £ is set 

here by the fluctuations of SU{i) Wilson lines in space. 
Perturbative corrections to that coefficient are small . 
Thus, we assume a Lorentz invariant structure The 
coupling ~ 22 between the chiral field and £ is chosen 
to reproduce m,r(T) This interaction term allows 

the Polyakov loop condensate to drive the chiral phase 
transition. While the inclusion of the chiral fields enables 
us to calculate familiar observables, it is otherwise largely 
passive. The behavior of i.e. confinement physics or its 
Z{2,) symmetry properties control the QCD transition. 
(In this regard, see also 0.) 

In our numerical simulations we solve the classi- 
cal Euler-Lagrange equations derived from the La- 
grangian (^) on a 64'^ space-like lattice. We employ a 
Robertson- Walker metric ds^ — dr^ — a? {T)A{TQrj)'^ — dx\ 
with scale factor a(r) = t/tq and a Hubble constant 
H = dlog(o)/dT = 1/r. Here, r and r] are the longi- 
tudinal proper time and the space-time rapidity, respec- 
tively We can vary the expansion rate by varying 
the initial time, tq. Further, we impose periodic bound- 
ary conditions and choose a lattice spacing ql = 0.25 fm. 
We draw random initial field values at each lattice site 
from a Gaussian distribution and then coarse grain the 
lattice afterwards to a correlation length of 2aL ■ The ini- 
tial temperature is chosen to be T(to) — 1. 01 Tc. For the 
chiral field fiuctuations, we assume |^ that {4>) = (</>) = 0, 
{(jp) = /A, and {(jP') = t;^/fm^ {v is the vacuum expec- 
tation value of the sigma field, a = (j)^). The i!-field is 
initialized in the deconfined minimum of the potential 
with real phase, {£) 0.7. The magnitude of the initial 
fluctuations, (£^) — {£)'^ w 0.1 was chosen such that the 
initial energy density is « 4Tc The results are es- 
sentially independent of the initial distribution of energy 
between fluctuations of £{x) and £{x). 

The presence of fluctuations modifies the potential in 
eq. (|l|). We renormahze it by expanding V{£, </)) at t = tq 
to second order in the fluctuations S£ and S(j) and then 
subtract the resulting SV{£,4>;S£,S(f>). This provides us 
with the renormalized equations of motion. 

We show the results of simulations for two differ- 
ent initial times; one appropriate for "slow" expansion 
(to = 2000 fm/c) and another {tq = 100 fm/c) for "fast" 
expansion. These two cases will suffice to illustrate the 
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radical changes which occur in the nature of the transi- 
tion as a function of the expansion rate. The slow expan- 
sion is, of course, many orders of magnitude faster than 
that describing the QCD transition in the early universe, 
while the fast expansion is slow on typical URHIC scales, 
and thus the indicators of non-equilibrium phenomena 
shown below should be even stronger in real URHIC. 




FIG. 1. \l\ field histograms obtained for the slow expan- 
sion. The times are a(r) = t/tq = 1.0443, 1.0456, 1.0467 and 
1.0480, from bottom to top in the left peak, respectively. 

The potential energy of the £-field in the co-moving 
frame is assumed to red-shift according to T*/Tq = 
l/a(T) = tq/t. However, the spectrum of long wave- 
length fluctuations on the lattice is determined by the 
dynamics of the phase transition; for example, in case of 
slow expansion with nuclcation and bubble growth those 
fluctuations on the lattice may not cool during the phase 
conversion. Only hard, short wavelength fluctuations 
which are integrated out and determine the coefficients 
in the effective potential (|l]), but are not treated explic- 
itly in the classical lattice simulation, cool according to 
the above prescribed T(r) law. 

Fig. |l| shows a sequence of histograms describing the 
distribution of \i\ for slow expansion. The transition is 
weakly first order and, accordingly, the effective potential 
for i exhibits two local minima for T in the immediate 
vicinity of Tc- Although {i{x)) is initially localized in the 
deconfined minimum at £q ^ 0.7, the fluctuations serve as 
a heat bath and trigger thermal activation over the bar- 
rier. Fig. shows that the confined and deconfined states 
coexist in different parts of the system. Well-defined bub- 
bles form and, as time progresses, grow, coalesce, and 
eventually complete the transition to confinement. 

Fig. H shows the corresponding histograms of \£\ for the 
fast expansion. Clearly, the double peak has disappeared, 
and the distributions merely broaden and move towards 
the confined minimum as the system cools. In this case, 
we conclude that the phase conversion happens through 
rapid decomposition. For a rough estimate of the relevant 
time scales, note that during a small time interval At, the 
hard modes cool hy\AT\/Tc ~ At/Atq. From classical 
nucleation theory a lower bound on the time scale 



for nucleation is T~^/'^ ^ l/Tc- For the model with 
b2{T) set by the string tension cr{T) supercooling by 
^ 1% suffices for the soft modes to "roll down" to the 
confined minimum Q rather than nucleation to occur. 
Setting AtTc ~ 1, we expect that for tq < 100 fm/c 
phase coexistence is not established. 




FIG. 2. \£\ field histograms obtained from the "fast" lattice 
simulation. The times are a(r) = t/tq — 1.0493, 1.0616, 
1.0739 and 1.0861 from right to left, respectively. 

Fig. H describes the time evolution of and its 

root mean square (RMS) fiuctuations. In the slow case, 
we see that the mean field makes the transition from 
the deconfined minimum to the confined minimum very 
smoothly. The corresponding RMS fluctuations increase 
immediately before phase conversion and illustrate the 
competition between the two minima. (For a second- 
order phase transition and very slow expansion those fluc- 
tuations would go "critical" [^). After the transition, the 
mean field shows no temporal oscillations and only small 
spatial fluctuations, indicating that the system is quite 
homogeneous after the transition. 
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FIG. 3. The mean field, {\£{x)\), (thick lines) and its RMS 
fluctuations in space (thin lines). The solid (dashed-dotted) 
curves describe slow (fast) expansion. 

By contrast, fast expansion leads to strong temporal 
oscillations of the mean field and large spatial fluctu- 
ations following the confinement transition. Thus, the 
scale of spatial homogeneity is smaller in this case. Also 
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note that the rapid phase conversion and decoherence 
seen for fast expansion lead to the violent and rapid de- 
cay of the Polyakov loop condensate into physical parti- 
cles. This, too, suggests explosive decomposition of the 
source. As discussed previously, strong event-by-event 
fluctuations of the pion mean-pt should emerge as a con- 
sequence Also, the "sudden" transition could help to 
understand the flavor composition in the final state . 

Fig. ^ summarizes the corresponding behavior of the 
chiral u field. The qualitative pattern here is the same, 
but the effects are enhanced. This indicates that substan- 
tial energy has been flushed quickly from the Polyakov 
loop condensate into the chiral field. This helps jus- 
tify our expectation of the approximate equality of out- 
ward and sideward HBT radii: The outward radius i?^ 
is larger than the sideward radius i?^ by an amount 
roughly equal to the variance of the pion production time 
(5r^ = (t^) — (t)^ ||lc| , pT] , [l8|| . For a smooth quasi-adiabatic 
transition, entropy conservation in a comoving volume cl- 
ement requires (5t^ to be large, while the non-equilibrium 
transition seen in our simulations is very fast. 
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FIG. 4. The mean field i(j(x)) (thick lines) and RMS fluc- 
tuations (thin lines) in MeV. The solid curves describe slow 
expansion. The dashed curves describe fast expansion. 

To summarize, real-time simulations of the Polyakov 
loop condensate model have demonstrated that radi- 
cally different phase transition dynamics can emerge in 
URHIC for realistically fast expansion rates when the 
effective potential changes rapidly in the vicinity of T^. 
The model considered is consistent with lattice QCD data 
for the pressure, energy density, and correlation length 
(or string tension) near T^- In real time, the picture 
that emerges is a confinement transition via spinodal de- 
composition which then triggers chiral symmetry break- 
ing. The timescales involved in URHIC seem much too 
short to admit the traditional near equilibrium transi- 
tion via bubble nucleation. As expected in any case of 
phase conversion via spinodal decomposition, we observe 
large RMS spatial fluctuations in the chiral fields after 
the transition, i.e. a small spatial homogeneity length. 
The near equivalence of outward and sideward HBT radii 
then follows as a natural consequence of the almost in- 



stantaneous decay of the Polyakov loop condensate into 
an "exploding" source of pions jlj] . Large event-by-event 
fluctuations of {ptj are another outcome of such a strong 
non-equilibrium transition. 
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FIG. 5. Additional figure: \l\ contours in the transverse FIG. 6. Additional figure: \l\ contours in the transverse 

plane, r) = Q, for the slow expansion, ro = 1970 fm. The plane, r) = for the fast expansion, tq = 98 fm. The times 

times are a(r) = r/ro = 1.045, 1.04875, 1.05031, from top to are a(r) = r/ro = 1.056, 1.069, 1.082, from top to bottom, 

bottom. Dark color corresponds to small I (confined phase) Dark color corresponds to small I (confined phase) while light 

while light color corresponds to large I (deconfined phase). color corresponds to large I (deconfined phase). 
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